Corrections to Tribimaximal Mixing from Nondegenerate Phases 
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Abstract 



We propose a seesaw scenario that possible corrections to the tribimaximal 

^ pattern of lepton mixing are due to the small phase splitting of the right- 

Q^' handed neutrino mass matrix, we show that the small deviations can be 

^1 expressed analytically in terms of two splitting parameters(5i and 82) in the 

D \ leading order. The solar mixing angle 612 favors a relatively smaller value 

' compared to zero order value (35.3°), and the Dirac type CP phase 5 chooses 
a nearly maximal one. The two Majorana type CP phases p and a turn out 

^ ' to be a nearly linear dependence. Also a normal hierarchy neutrino mass 

O , spectrum is favored due to the stability of perturbation calculations. 
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I. INTRODUCTION 



Due to the out-standing achievements in the neutrino oscillation experiments [1-18], we 
have convincing evidences on neutrino mass and lepton mixing. In the weak basis that the 
charged-lepton mass matrix is diagonal, real and positive , we assume massive neutrinos to 
be Majorana particles and parametrize [19,20] the lepton mixing V as 



AS « « c. c. c. ^*<5 



"512^23 <^12'^23"5l3^* C;^2C23 '^12'^23'5l3^* '^23'^13 6*°" . (1-1) 



V "5l2'^23 ''-^12'^23'^13^"' ^12^23 "5l2C23"5l3^*'' <^23'^13 




Now the global fit [20-22] of current experimental dada yields sin^ 6'i2 ~ 0.304lo:oi6 > 
sin^ 023 ~ O.SOloloe )Sin^ ^13 < 0.016 at 1 a. So the so-called tribimaximal mixing [23-26] 
pattern is an excellent approximation for these physical values, which appears as the form 
of 

/ 2/V6 1/VS \ 
= 1/V3 l/x/2 . (1.2) 

V 1/V6 1/V2J 

It corresponds to sin^ 612 = 1/3 , sin^ 623 = 1/2 and sin^ ^^13 = in standard parametrization. 
This simple form with small integers motivates neutrino theorists to consider some under- 
lying structures. The most intriguing one is flavor symmetry among generations, especially 
for discrete groups [27,28] such as A4 , ^'3 , Z2, which can give predictable values or/and 
relations of mixing parameters. In general, nearly tribimaximal mixing is also allowed, so 
small corrections to the standard form in (1.2) arc interesting and necessary. Many sources 
can give these corrections such as charged lepton sector contributions [29-31] and renormal- 
ization group effects [32-34] . 

On the other way, seesaw mechanism [35-39] can naturally explain the smallness of three 
left-handed neutrino masses. In the simplest type I framework, the effective Majorana mass 
matrix of neutrinos Mj, is related to the Dirac mass matrix Md and the heavy right-handed 
Majorana mass matrix Mr by the relation 

M^c^MdM^MI. (1.3) 

To explore the structure of the seesaw formula, the effective mixing can generally be derived 
from both Md and Mr. One of them can be chosen to be diagonal or even identity matrix 



for simplifications. One kind of models in ref. [40-42] takes a diagonal Md for instance. 
Meanwhile, another kind of seesaw models [43,44] starts from a unit form of Mr. 
Now for the small corrections of tribimaximal mixing, another source [43,44] in realistic 
seesaw models attracts our attentions. The leading order mixing matrix is totally derived 
from the Dirac mass matrix M^, and Mr has the simplest form of MqI in the symmetry 
limit (Mq being a common mass scale and / the identity matrix) . The departure of Mr 
from the unit matrix gives small corrections to the effective mixing matrix. There are two 
methods in Mr to break the unit form, one is the nondegenerate masses, another is the un- 
degenerate phases: the phase breaking. In this letter, we want to discuss the second case, we 
consider small phase splitting for complex matrix Mr, and maintain the degeneracy of the 
right-handed neutrino masses. We find that besides small corrections of the three mixing 
angles, we get nontrivial CP violation phases. 

The remaining part of this paper is organized as follows. In section II, we talk about the 
realistic models from phase breaking, giving the predicted values of mixing parameters. In 
section III, we do some numerical analysis, displaying the correlations of these parameters. 
We conclude our topic and give some remarks in section IV. In the Appendix, we give our 
calculations from perturbation approximations. 



In our model, we adopt the basis started in [43,44]. we think M^ can be diagonalized by 
tribimaximal mixing matrix in the manner of VQ~M£,VQ=Diag{x, y, z} which is constrained 
by a discrete flavor symmetry(such as symmetry), (x, y, z) are positive mass parameters. 
For Mr, we take the form of MqI in the symmetry limit, obviously it is invariant under the 
symmetry. Then we have 



In the following, we want to give small phase splitting to Mr, and derive the light neutrino 
mass spectra and small corrections for lepton mixing parameters compared with the leading 



II. REALISTIC MODELS 



/x 

M,^MdM^^MI=Vq y 

yo z 




(2.1) 



3 



values of Vq including CP violating phases . 



A. Scenario A 



In this scenario, we give Mjn i a small phase compared to the others, then Mr take the 
form of 



-iSi 



Mr = Mo 

Similar to (2.1), the effective neutrino mass matrix Mj, turns out to be 





V 




(2.2) 
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M,^MdMr'mI=Vo 
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2 / (3 + 2ei)uj^7f V2ujri^ 



3Mq 



0^ 



V 







(3 + 61)77^ |yo^ = yoMWVo^, 
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(2.3) 



with u = x/y , 1] = y/z and ei = e^^^ — \ . The corrections to Vq come from the diagonalization 
of matrix M^^^. Assuming M^^ = V^M^^V^ , where M^, = Diag{mi,m2,m^} with mj(i = 
1,2,3) being the Majorana neutrino masses. When we parameterize as 



/ cos a sin ae 
- sin cte*^ 

V 








\ 




(2.4) 



we can derive these mixing parameters and masses by the relation of M^ = V^M^^^V{ 
in terms of ei , ui , f] , z and M^ . In our calculation, we view ei (or equivalently 5i) as a 
small parameter, and give results with the leading terms of the power series of Si , such as 
ei = e^^^ — 1 ~ —Sl/2 + i5i . By the perturbation calculations, we derive the relation of 
cos/3 ~ l/6sin/3(5i , immediately we have cot /3 ~ 5i/6 . Other parameters are listed below 

2^/2 



tan 2a ~ 



3 l+uj- 



(2.5) 



mi 



z'^uj'^rf 



Z^Uj'^Tf 



Mo ' 9 ' 91 + LJ^ Mo 



1712 



- Mo^ 9^'^ 91 + 00^^'^- Mo^ g'^^ 



(2.6) 
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and 



2 1 
m3 = — — , tan2(7i ~ -8\ , tan 2(72 — 77(^1 , tan 2(73 
Mn 3 3 



0. 



(2.7) 



For the second expressions of m\ and , we omit the contributions of higher powers of 
Lo . Combining together with Vq and making a rephasing transformation, we can get the 
standard expression of neutrino mixing matrix V as eq.(l.l) 




coso; sinoie 0\ 
'■^ cos a 

ly 



sm ae 




/ e'"^ 


V 



\ 








r,l<73 



12 



/ COS 9 

sin ^12 



12 



sin^ 
^ cos ^12 



\ /e^p 
1 



V2 





V 



0\ 

e'^ 



(2.8) 



Ij 



V 72^i"^i2 -73Cosei2 

where pi{i = 1,2,3) can be rotated away by redefining the phases of three charged-lepton 
fields. Then the standard parameters appear as 

I r , a; , 1 



sin 9 



12 



:{1 



and two standard Majorana phases are 



(2.9) 



2,^ 2u 
tan2p~ -{1- )Si, 
3 1 + o;"^ 

1, 2u; 
tan 2(7 ~ -fl - ^ -)5i . 



(2.10) 



3' l + u"^' 

In this scenario, we get a small deviation from tribimaximal mixing for 9i2 , and two Majo- 
rana phases, but ^13 and ^23 remain unchanged. Because of the vanishing value of ^13 , we 
cannot obtain the information of Dirac CP phase. Contrarily if we give Mij2,2 another small 
phase compared to Mij3 3, we can obtain a non- vanishing ^13 and then a non-trivial Dirac 
CP violation phase. 



B. Scenario B 



As discussed in the end of last section, we take a general form of Mr as 

-iSi ^ 
Mr = Mo\ 
1 



(2.11) 
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Then we will repeat the same procedure as before. Firstly the effective neutrino mass matrix 
are 

fx 0\ . /e^^i 0\ fx 0\ 
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^2 / (6 + 4ei + e2)cj V (2ei - e2)y2cu?72 -^^2^,,^^ 

(2ei - e2)V2u7f (6 + 2ei + 2e2)rf VQe2V | V^j"' = V^M^^^V^ . (2.12) 
-A/3e2CJ?7 A/6e2r7 6 + 3e2 



6Mo 



Then we define M^^ = VgM^Vg^ , with V2 unitary matrix parameterized as 



^ Di as 0L2 ^ 
\-ai, -al D^j 



\ 



(2.13) 



V e'"^ ) 

where aj = sin ^^e"*'^-' (for j — 1,2,3) , Di — ^^1 — |Q;2p — IctsP and similar definitions for 
£)2 and £)3. By using = V2^M^^)V2* , we arrive at the approximate expressions for 
6i and 0j via perturbation calculations, (we will give a careful calculations of perturbation 
approximations in the Appendix.) 



sin 6*3 sin 03 ~ a;(2(5i — ^2) ,sin6'2sin02 — — 'jjr\^2 , sin 6*1 sin 0i ~ 77^2 
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(2.14) 



and 



\/2 

sin 6I3 cos 03 ~ - (5i52 + (^2) , 

io 

\/3 

sin Q2 cos 02 — rr^- <^ ^ (^2(2(5i - (^2) , 
lo 

sin cos 01 ~ -^^("^^ - 77^)^2(251 - ^2) • 



(2.15) 



And the masses are 



2:^0; V /. 1 c-2 1 c- c- 5 ^2n 



777,2 



'0 
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18 



Finally we give the phase-angles Ui in eq.(2.13) 

tan2(7i ~ ^(45i + ^2) , tan 2(72 — ^((^1 + (^2) , tan 2(73 — ^<^2 



(2.16) 



(2.17) 
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A normal hierarchy mass spectrum(small values of cu and 77) is favored in this scenario due to 
the stability of perturbation calculation (see the Appendix). So for the analytical expressions 
in (2.14) (2.15) and (2.16), we express them with the leading terms of the power series of u; 
and f]. More detailed expressions and discussions can be found in the Appendix. 
Now we want to exhibit the mixing matrix V — V0V2 with the standard parametrization as 
in eq.(l.l). By a rephasing transformation, we get 

\ 



V 

( 



1 1 1 

, J- L A. I 



V6 

i_ 

'6 



-at 



"3 

D2 



i5 



^12^2Z ^12*23*13^ 
"•^12*23 ~ '^12'^23*13^ 















*23'^13 



(2.18) 



V ^12*23 ~ C22C23S]^3e'" "^^2*23 ~ '^12'^23*13^'" '^23'^13 / 

By a lengthy but straightforward calculation, we can get the standard predictions. The 
Dirac type CP phase is predicted as 



r 5 uj 

cos () ~ 

61 — 



UJ 



■(25l-52), 



(2.19) 



and the Majorana CP phases are 



tan2p~ -{l-2uj){25i- 82) , 
3 



tan 2cr ~ -(1 
6^ 



2uj){25i-52). 



(2.20) 



About the mixing angles, we can obtain a small departure from the tribimaximal mixing in 

the manner of 



sin 6 



13 



6 ^l + ry2 



curj 



1 + u'^rf 



tanei2 ^ ^{1 - - S,52 + Si)} , 

2 6 



tan^ 



23 



■u;r]52{2Si - 62) 



(2.21) 



Immediately, when ^2 ^ 0, the predicted values return to the corresponding ones in Scenario 
A as we expect to, and the Dirac phase 5 is undetermined due to zero of sin ^13. The small 
parameter 82 is responsible for the non- vanishing ^13 and non-maximality of ^23, and the 
mixing angle 612 depends on both 5i and 52- In this scenario, we predict a nearly maximal 



Dirac phase, so the rephasing invariant Jarlskog parameter J [45,46] is determined only by 
the three mixing angles in leading order 

J ~ ^ sin ^13 - ^?7(1 - ^)l<^2| . (2.22) 

Indeed, the most general Dirac mass matrix can be diagonalized by two distinct unitary 
matrices Vl and Vr, so if we identify Vl as the exact tribimaximal form, the high order 
corrections come from the contribution of Nr = V^M^^V^, which is a unitary symmetric 
matrix if the mass parameters in Mr are degenerated. This case has been discussed for the 
(light) effective Majorona mass matrix with exact degenerate masses [47,48]. They found 
there are only two mixing angles and one CP phase in the parametrization of mixing matrix, 
but it does not include the case with the same CP parity as a limit. In our scenarios, we want 
to get higher order corrections from the diagonalization of , so we are only interested 
in the nearly unit form of Nji . It is convenient to parameterized as = V^Ds- where 
Vjv is a real unitary (orthogonal) matrix and D^^ is a diagonal unitary matrix with small 
denoting the small departure from unit form of A^^. In the beginning of our calculation, we 
have identified V^v with the explicit form of Vq for simplicity. Although this is only a special 
example, but it is enough to reveal the property of our scenarios and correlations of these 
corrections. In next section, we will show them numerically. 

III. NUMERICAL ANALYSIS 

Scenario B returns to Scenario A when ^2 ^ , so we can only concern the general case. 
To do numerical analysis, we will include higher order contributions of u and r^as in ((5.2), 
(5. 7) -(5. 9)). Firstly, the correlations of cu and rj come form the ratio of Aml^; and Am%^ 
by the relation of 

Using the numerical results in [21,22] which reveal that Aml^; = 7.65lo;|o x lQ~^eV^ and 
Am\f^ — 2.40lo!ii x lO~^eV'^ at 1 cr , we can obtain the allowed region of u; versus rj in 
FIG.l. The mass spectrum changes from hierarchy region to nearly degenerate region as 
cu grows. And rj takes a nearly fixed value in the hierarchy region. As discussed in the 

8 



Appendix, only the normal hierarchy mass spectrum is valid in our scenarios, so we restrict 
the parameter rj < 0.5 in our analysis. 

The departure of mixing matrix from the standard tribimaximal mixing depends on the 
two small parameters Si and 62- If we change the values of 6i , we can get the correlation 
relations between mixing parameters (mixing angles and CP phases) defined in (1.1). Scan- 
ning Si within a reasonable range( [-0.5 , 0.5 ]), we give five pictures FIG.2-FIG.6 as typical 
examples. The angles and phases are measured in degrees. Some comments are listed below. 

• For the relation of between 612 and ^13 in FIG. 2 , we can see that the correction of 612 is 
only at the left side, which has been revealed in (2.21) for leading corrections, indicating 
that a relatively small 9i2 is favored, which is in accordance with the best fit point (33.2°) 
of global analysis [?]. The magnitude of ^13 is too small (< 2.5°) to be measured in the near 
future. The sensitivity of the proposed reactor neutrino experiments to ^13 is at the level of 
^13 ~ 3°(sin2 2^13 ~ 0.01) [20,57]. Moreover, a larger deviation of 612 from 35.3° (sin"^ ^J) 
implies a more stringent constraint on ^^13, including a lower bound. 

• In FIG. 3 , the deviation of 823 can extend to both side, but with an unsymmetrical shape. 
The magnitude of this deviation is less than 1°. A larger ^13 allows a smaller region of ^23- 
Oppositely, a larger deviation of 623 from maximality implies a smaller range for ^13. 

• The Dirac type CP phase 5 favors a nearly maximal value which can be understood in 
(2.19) and in FIG. 4. The relatively broad width results from the large range of {26i — 62)- 
The figure shows that the more deviation of S from maximality, the more stringent upper 
bound ^13 suffers. The Jarlskog parameter J can reach the order of 10~^ at most, which is 
hmited by the relatively small value of ^13 . 

• Just as revealed in FIG. 5, the two Majorana CP phases p and a have distributions among 
[—180°, 180°] peaked at zero, which is consistent with leading predictions of (2.20). Also, 
FIG. 6 shows that p and a are strongly correlated with each other, having a nearly linear 
dependence. The slope of p versus a defined in (2.20) approximates to 0.5 , and the corre- 
sponding one in FIG. 6 hes between 0.75 ~ 1.0, where the higher order contributions of u; 
and 77 have been included. 
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IV. CONCLUSIONS 



To conclude, we discuss possible corrections to the tribimaximal pattern of lepton mixing 
in a simple seesaw model, which is the results of phase breaking of right-handed neutrino 
matrix. We consider small phase splitting for complex matrix Mr, but maintain the de- 
generacy of the right-handed neutrino masses. The breaking of the unit form of Mr in 
the seesaw model gives small corrections to the zero order form of the mixing matrix. As 
revealed in (2.21), the corrections of ^13 , 9i2 , and ^23 are of order 0{5i) , 0{5l) , and 
respectively. 9 12 favors a smaller value compared to the zero order value of 35.3° . The Dirac 
type CP phase 5 is hkely to have a nearly maximal value. For the two Majorana CP phases 
p and (T, they have a nearly linear correlation. 

A normal hierarchy mass spectrum is favored in our scenarios, leading to insignificant renor- 
mahzation group effects [32-34,43,44] of the mixing matrix. The exact mass degeneracy of 
the righthanded neutrinos forbids CP violation in the lepton-number- violating decays [49], 
so there is no thermal leptogenesis [50]. Generally, we should include both effects men- 
tioned in the introduction: the nondegenerate masses and the nondegeneratc phases. But 
our discussions are valid when the mass degeneracy breaking is much smaller than the latter 
one. The general cases including both (and the corresponding leptogenesis) are certainly 
interesting and need further discussions. 

There are other interesting aspects on the corrections to tribimaximal mixing based on either 
different assumptions [51,52] or general parameterizations [53-55]. They can only be distin- 
guished by precision measurements in future reactor and accelerator neutrino experiments 
[20,56,57]. 
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V. APPENDIX: PERTURBATION CALCULATIONS FOR DIAGONALIZING 

THE NEUTRINO MASS MATRIX 



When we want to diagonalize the matrix of M^^^ in (2.12), we refer to the assumption 
that parameters Si are small enough to do perturbation calculations. In the relation of 
— V2'Mj;^^^V2 ,we think equations are reahzed in each order for the power series of 
5i. Then when we put the approximation to the order of 6i {0(61,62)), v/e can get three 
equations by the relations of three off-diagonal elements 

iV2{25i — 62)007^ + sin ^3(cos 03 + i sin (f)3)6u}'^r]^ — sin ^3(cos (ps — i sin ^3)677^ = , 
—i\/?>620Jri + sin6'2(cos02 + ism(i)2)Quj'^'rf — sin6'2(cos02 — isin02)6 = , 

iVQS2r) + sin ^1 (cos 01 + isin0i)6?7^ — sin^i(cos0i — isin0i)6 = . (5.1) 

Then we get the results as 

sin ^3 sin 03 ~ - „ .^^^ (2(5i - ^2) , 
6(1 + cu^j 

sm^2sm02 ^ ^7— — ^^^2 , 
6(1 + uj^r]'^) 

a/6 

sin6lisin0i ~ -— — — -62, (5.2) 
6(1 + 7]^) 



and 



sin 63 cos 03 = sin 62 cos 02 = sin 62 cos 02 = . (5.3) 



Which indicate maximality for 0i(cos0j — 0). We want to know the departure from max- 
imality in order of 0{5i,52), so terms in order of 0{Sf ,6162,62) must be considered in 
equations of (5.1). Up to this order those equations have the form of 

iV2{26i — 62)^07}"^ + i sin 63 sin 03677^(1 + o;^) + sin 63 cos (p^Grj^luj^ — 1) 
+V2{-6f + -j-)urf - sme^sin (t>3rf[{A6i + 62)^'' + (25i + 262)] 
—6 sin 9i sin 0i sin O2 sin 02 — \/6 sin 62 sin 0252^7 + \/3 sin 61 sin 4>i620Jrj = , (5.4) 

—i\/362UJr] + i sin 62 sin 026(1 + lj'^i]'^) + sin 62 cos 026(0;^?]^ ~ 1) H — tt^I^V 
— sin Q2 sin 02[(45i + 62)oJ^rf' + 352] — 677^ sin 9i sin 0i sin 6^ sin 03 

— -\/2 sin^i sin0i(25i — 62)007]"^ — v^sin 63 sin 0352?7 = , (5.5) 

11 



\/6 

i\/6S2rj + i sin 6i sin 0i6(l + r^^) + sin 62 cos 026(?7^ ^1) 2^^^^ 

— sin ^1 sin 0i[(25i + 262)1]^ + ^62] — QuJ^rf' sin 62 sin ^2 sin ^3 sin ^3 

— -\/2 sin^2 sin 02(2^1 — 82)0^7^ + -\/3sin^3 sin (f)^520or] — . (5-6) 

Firstly, the imaginary parts of the equations are in order of 0{Si, S2), giving the same values 
of sin^jsin^i as in (5.2). Contrarily, from the real parts of the order of 0{5j,6i52,52), we 
can obtain the expressions for sin 9i cos 0j . These results read as 

sin ^3 cos 03 - - 3g^^'^^4^ {2(l-a;^)^i^ - 2{1 - u;')S,S2 

-(1 + 2(^)^2 + 3 — /I , 2^n ^ 2 2N — ^2}, (5.7) 

sin 6*2 cos 02 — — V^*^^ ^^ 62(261 — (52)|— cu^r?^ 
^ 18(1 -cjV) 

, (l + a;V)(l + r7^ + cuV) ^ 

+ (l + u;2)(l + r^2) i' ^^■^) 

sin 6*1 cos 01 ~ -— — TT(^2(25i - 52){-77^ 

36(1 — Tj^) 

.2 (l+^^)(l+^^ + ^V) l 

(l + c.2)(l + ..V) ^ ■ ^ ^ ^ 

Following the spirits of perturbation calculation, higher order contributions only give small 

corrections, meaning that the results are stable. Comparing these results with those of 

(5.3), we can see that if the parameter uj or 77(or both) is very close to 1 (quasi-degenerate 

or inverted hierarchy mass spectrum), the property of perturbation will be ruined. So a 

normal hierarchy(small uj and rj) neutrino mass spectrum is favored in this scenario. Then 

for analytical expressions, we can simplify them with the leading terms for power series of 

UJ and T] . For the results of sin 9i sin 0j and sin 9i cos 0j , we have 

sin6'3sin03 ~ -^uj{25i -82) , sin 6*2 sin 02 — -^'^^(^2 ,sin6'isin0i ~ -^r)S2 ; (5.10) 



and 



V2 

sin 9s cos 03 ~ -—uj{Sl - 6162 + Sj) , 

io 

sin 6*2 COS 02 ~ —ujr]62{26i - ^2) , 
lo 

\/6 

sin^i cos 01 ~ -i7rry(a;' - ^^2(261 - S2) . (5.11) 
3d 
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Similarly, for the diagonal elements in M^^ — Vg^M^^^V^* , we can get the information of rrii 
and CTj as 

zW 1 1 + 3c^2 1 1 + 3cu% ^ 1 1 + 4:7]^ + 3loY ,2^ 

-3 - ^(1 - i^.^ + 1,. 2+.y3.V _ 
Mo^ 8^ 12 ' (l + 772)(l + a;V) 

and 

tan2cTi ~ ^(45i + ^2) , tan2(T2 ^ ^(^i + ^2) , tan2(T3 ~ ■ (5.13) 

D O Zi 

For small mass ratios, the three masses are expressed as 

z'' 1 



^3^xr(i-o<^2)- (5-14) 



Mo 
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FIGURES 




FIG. 1. Allowed parameter space of lo and r] from la region of the two Am^ 
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FIG. 2. Allowed parameter distribution of 612 and 613 for {61 62) within [-0.5,0.5] 
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FIG. 3. Allowed parameter distribution of 6*23 and ^13 for {61 82) within [-0.5,0.5] 
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FIG. 4. Allowed parameter distribution of 6 and ^13 for (61 62) within [-0.5,0.5] 
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FIG. 5. Allowed parameter distribution of p and for {5i 82) within [-0.5 ,0.5] 
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FIG. 6. Parameter dependence of p and a for {5i 82) within [-0.5,0.5] . 
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